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Abstract: Non-equilibrium black hole horizons are considered in scaling theories 
with generic Lifshitz invariance and an unbroken U(l) symmetry. There is also a spe¬ 
cial form of hyperscaling violation associated with a non-trivial conduction exponent. 
The boundary stress tensor is computed and renormalized and the associated hydro- 
dynamic equations derived. Upon a non-trivial redefinition of boundary sources 
associated with the U(l) gauge held, the equations are mapped to the standard 
non-relativistic hydrodynamics equations coupled to a mass current and an external 
Newton potential in accordance with the general theory of |Q. The sheer viscosity 
to entropy ratio is the same as in the relativistic case. 
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1. Introduction 


The AdS/CFT correspondence ^ ^ relates the anti-de Sitter (AdS) space-time 
to conformal held theory (CFT) on the boundary. It gives a semiclassical description 
of strong coupling physics in the dual held theory in terms of string theory or its 
low-energy limit: (super) gravity. At hnite temperature and in the long wavelength 
regime, the dual held theory can be ehectively described by huid mechanics and it 
can be related to black holes in AdS space-time. The huid/gravity correspondence 
was hrst studied using linear response theory i, I, 0, i- 

Subsequently, fully dynamic descriptions were studied using boosted black holes 
in asymptotically AdS geometries that led to relativistic huid dynamics in the dual 
CFT In this formalism, the huid variables are encoded in the near-equilibrium 
black hole solution and the huid equations appear as constraints on the solution 
imposed by the bulk equations of motion. 

Recently, generalizations of the AdS/CFT correspondence to theories with non- 
relativistic scaling symmetry have been studied. In particular, many condensed 
matter systems have critical points with non-relativistic scale invariance m 0 
Some of these systems have Lifshitz or Schrodinger symmetry 


12, 13, 14 


17|. Moreover, the hydrodynamics of charge and energy in such systems may be 


interesting as has been argued recently for the case of cold fermions at unitarity, 
p!8[| , other strongly correlated systems, []^ and graphene, 0. 

Holographic techniques have been generalized to geometries with Lifshitz or 
Schrodinger symmetry in connection with applications to condensed matter systems 
[PT| , P2| , |25| , ^ In particular, in [^, ^ |^, all quantum critical holo¬ 

graphic scaling theories with a U{1) symmetry respecting translation invariance and 
spatial rotation invariance were classihed in terms of three scaling exponents. Two 
of them {z, 9) appear in the metric while another exponent, which is referred to as C, 
2^ ^ 1^, appears in the prohle of the U{1) gauge held.^ The exponent z is the 


m 


Lifshitz (dynamical) scaling exponent, and 6 is the hyperscaling violation exponent, 
mM- Even though such theories have been studied intensively many of the aspects 
are still unclear and, in particular, hydrodynamics with Lifshitz scaling symmetry is 
not fully understood. 


related work was presented in Q. 

^This charge exponent control the anomalous scaling of the charge density, even if it is conserved, 
has also been introduced independently in Q and was studied in more detail in and [^). 
The reason for the existence of anomalous charge exponent despite conservation is the RG running 
of the bulk coupling for charged degrees of freedom. 
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More recently, it was found that the boundary theory dual to space-times with 
Lifshitz asymptotics can be described in terms of the torsional Newton-Cartan grav¬ 
ity theory, which is a novel extension of the Newton-Cartan gravity with a specihc 
torsion tensor. The application of the Newton-Cartan theory to non-relativistic con¬ 
densed matter systems (namely the Quantum Hall effect) was first discussed in . 
Interactions between the torsional Newton-Cartan gravity and matter were discussed 
in [CT|. The correspondence between the Lifshitz space-time and boundary torsional 

for specihc Lifshitz geometry and 


Newton-Cartan theory was hrst found in pl3| , 
further studied in ISSl, [ 


37, 38]. In these works, the correspondence is studied by 


using the vielbein formalism, in which an appropriate combination of the vielbeins 
and bulk gauge helds is considered. It turns out to be very useful to use vielbeins to 
study the boundary theory. This is consistent to the holographic renormalization in 
the asymptotically Lifshitz space-time, in which the scaling dimension is calculated 
by using the vielbein |^, Q. Counter terms in Lifshitz space-time were discussed 
in generality in by using the Hamilton-Jacobi formalism. 

The huid/gravity correspondence for non-relativistic huids has been studied in 
^ for a special case of the Schrodinger geometry which is related to ordinary 
AdS by the TsT transformation. In these studies, the non-relativistic fluids are 
obtained by the light-cone reduction of the relativistic fluids. The generalization to 
the charged fluid case is studied in . 


For the Lifshitz space-time, the correspondence to relativistic fluids with Lif¬ 
shitz scaling, in which the velocity field is dehned by a normalized lorentz vector, 
was studied in P5| , ^ |T7| . In these works the anisotropic direction of the Lifshitz 
symmetry depends on the frame. The fluid appears on the surface at finite radius or 
on the horizon, contrary to the Newton-Cartan theory which appears on the bound¬ 
ary. The hydrodynamics found contains an antisymmetric part in the hydrodynamic 
stress tensor that contributes a new transport coefficient to the dynamics. 

In this paper, we consider the fluid/gravity correspondence for Lifshitz geome¬ 
tries and the relation to fluids in boundary non-relativistic theories with Newton- 
Cartan symmetry. We consider black holes in Lifshitz space-time with unbroken 
U{1) gauge symmetry that are solutions of the Einstein-Maxwell-dilaton (EMD) 
theories. Although, the geometry has Lifshitz scaling symmetry with dynamical ex¬ 
ponent z, the bulk solution has hyperscaling violation^ due to a nontrivial exponent 
conduction exponent associated with the gauge field and the non-trivial running of 
the dilaton. 

We consider the black-hole solution of the theory, boost it using Galilean boosts 
and then we make all parameters of the solution including the velocities, x-dependent. 
We then proceed with the standard analysis introduced in 0: we solve the bulk equa- 


^This is distinct from what is called hyperscaling violation in condensed matter physics. Our 
definition of generalized hyperscaling violation is based on the existence of scaling but also the 
existence of anomalous dimensions in the charged sector. 
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tions of motion order by order in boundary derivatives and compute and renormalize 
the (fluid) stress-energy tensor. We also calculate the entropy current and consider 
the thermodynamic relations. What we hnd is as follows: 

• The stress-energy tensor is expressed in terms of the fluid variables: velocity 
held u*, energy density £ and pressure P, but also contains the (particle num¬ 
ber) density n and external source Ai associated to the P(l) symmetry current. 
It satishes the condition for Lifshitz invariant theories z£ = {d — \')P. 

• By comparing the stress-energy tensor and the constraints in the bulk equations 
of motion, we hnd that the conservation law of the stress-energy tensor is diher- 
ent from that of relativistic theories but agrees with that in the Newton-Cartan 
theory. In the stress-energy tensor required a modihcation (improvement) 
in order to satisfy the trace Ward identity. Our stress-energy tensor satishes 
the Ward identity without such a modihcation. 

• The role of the (unbroken) U(l) symmetry in this class of theories is important. 
It should be noted that this U(l) symmetry is responsible for the Lifshitz 
background bulk solution. We hnd that it behaves very closely to the U(l) 
mass conservation symmetry in non-relativistic hydrodynamics. 


The huid here is non-relativistic and this is diherent from the relativistic huids 
analysis in [^, Even though the continuity equation and energy conser¬ 


vation equation agree with those in the ordinary non-relativistic huids, the 
Navier-Stokes equation is diherent from that in the ordinary non-relativistic 
huids. The ehects of pressure become much larger than other contributions 
and some terms with velocity held in the Navier-Stokes equation are absent 
in our result. These absent terms are replaced by external source-dependent 
terms associated to the U(l) gauge held. 

By redehning variables and allowing a (Milne-invariant) Newton potential in 


our sources, 


we can map the huid equations to the standard non- 


relativistic huid equations coupled to the torsional Newton-Cartan geometry 
in the presence of a Newton potential. This is a universal result that we hnd 
interesting and far-reaching 


Even if the gauge held is not identihed to that in the Newton-Cartan theory, our 
huid can be interpreted as a non-relativistic limit of a huid with Lifshitz scaling 
symmetry. In ordinary non-relativistic huids the pressure is much smaller than 
energy density but in our case they are at the same order due to the Lifshitz 
scaling symmetry. For non-relativistic huids with Schrodinger symmetry the 
huid equations are obtained by introducing the light-cone dimensional reduc- 

Instead, our non-relativistic limit arises naturally as a 


tion in 


rather ordinary limit. 
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• We find that the form of the fluid equations is independent of the Lifshitz ex¬ 
ponent z as well as of the (non-trivial) conduction exponent. It is only the 
constitutive relations (equation of state) that depend on these scaling expo¬ 
nents. 

• The entropy satishes the local thermodynamic relation with the energy density 
and pressure. The divergence of the entropy current is non-negative, compatible 
with the second law. 

This paper is organized as follows. In Section ||, we introduce the model and 
its solution of Lifshitz space-time. Then, we hrst focus on the case with scaling 
exponent z = ‘1. In Section we introduce the hydrodynamic ansatz. In Section ^ 
we solve the equations of motion by using the derivative expansion and obtain the 
solution to the hrst order. In Section |^, we calculate the stress-energy tensor on the 
boundary and study its conservation law. In Section we investigate the entropy 
and thermodynamic relation. In Section H, we consider the generalization to general 
z. Section |] is devoted to conclusion and discussions. 

Appendix ^ contains a list of the variables used in this paper and their dehnition. 
It contains also comparison of variables with two other papers in the literature. More 
details on the calculation for hrst order solution and boundary stress-energy tensor 
are described in Appendix P and [p, respectively. In Appendix 0, we consider the 
analysis of the general counter terms. In Appendix |^, we discuss the regularity of 
the gauge held at the horizon. More details on the solution and counter terms for 
general z are discussed in Appendices ^ and B respectively. 

2. U(l)-invariant, hyperscaling-violating Lifshitz theory 

We consider a holographic theory with Lifshitz scaling and an unbroken U(l) global 
symmetry. The dual {d -|- l)-dimensional gravity theory will have a massless U{1) 
gauge held and a dilaton 0. The bulk action is given by 

s = (^R-2A- . (2.1) 

where F = dA is the held strength of the gauge held, and A is a dimensionless 
coupling constant of the bulk theory. The equations of motion are given by 

Rfxu = {^F^pF/ - ( 2 . 2 ) 

0 = V^(e^‘^F^'^), (2.3) 

□0 = . (2.4) 
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This model has the Lifshitz geometry as a solution; 

ds^ = H- Y + r‘^{dx'Y, 

i 

with the following the gauge field and dilaton; 

A, = , eA* = , 


(2.5) 


( 2 . 6 ) 


where the parameters z, a and /i are related to the parameters of the action (coupling 
constants) as 


A 


2 


d — 1 
z — 1 ’ 

{z — d — l){z + d — 2) 


2 


(2.7) 

(2.8) 


fior — 


2{z-l) 


z d — 1 

This solution, although well known in the context of cosmology 
was first studied in holography in ^ and was generalized in 
The metric (|2.5|) has the Lifshitz scaling symmetry 


(2.9) 

since a long time, 


t ^ c^t , X* —)■ ex* , r ^ c , (2.10) 


but charge density is transforming anomalously, violating hyperscaling. The reason 
is the running of the dilaton. We define hyperscaling violation here in more general 
terms as the presence of scaling but the failure for at least one physical observable to 
scale with standard dimension d scaling. The solution above is a solution with hy¬ 
perscaling violation, [^, ^|, although the hyperscaling violation exponent 6 coming 
from the metric vanishes. 


This model also has a black hole geometry as a solution p4 


where 


ds^ = -r^-f{r)dt^ + 


/ = 1 


f{r)r^ 


z+d-l 
_0 _ 

Y^z~\~d — 1 


^r^(dx*)^ 


The Hawking temperature of the black hole is given by 


T = 


An 


( 2 . 11 ) 


( 2 . 12 ) 


(2.13) 


The gauge held and dilaton take almost the same form as in the the zero temperature 
solution 


At = a(r"+'^-i - r^+'^-^), (2.14) 

but At must vanish at the horizon to be regular. 

In this paper, we focus on the case of d = 4. Extensions to other dimensions are 
expected to be straightforward. 
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3. Hydrodynamic ansatz 


In this section, we introdnce an ansatz for the geometry which describes the physics 
of fluids in the boundary quantum held theory. We use the method proposed in p. 
We will also hx the Lifshitz exponent to be z = 2. In a later section, we will discuss 
other values of z. 

In order for the regularity at the horizon to become evident, we change to 
Eddington-Finkelstein coordinates; 

+ 2rdt+dr + r‘^{dx^Y ■> (3-1) 

where the null coordinate is dehned by 

dr 


dt+ = dt + 


r^f 


The gauge held becomes 


A = a[r^ — Tq) — ar‘^dr , 


(3,2) 


(3,3) 


where we have hxed the Ar = 0 gauge in the Feherman-Graham coordinates. Here¬ 
after, we always take the Fddington-Finkelstein coordinates and t will stand for the 
null coordinate 

To implement the hydrodynamic ansatz, we hrst boost the black hole geometry. 
In the case of the ordinary Schwarzschild-AdSs, the boundary held theory is a rela¬ 
tivistic conformal held theory, and hence the Lorentz boost is employed. The Lifshitz 
geometry, however, corresponds to the (torsional)-Newton-Cartan theory, [^ -[08 


Therefore, we perform a Galilean boost on the black hole geometry. The metric 
becomes 


= —{r'^f — v‘^r‘^)dt^ -|- 2rdtdr — 2r^n*dt d;r* -f- r'^^dx'Y 
The gauge held and dilaton are not ahected by the Galilean boost; 


(3.4) 


A = a{r^ — Tq) dt — ar'^dr , e'^'^ = pr (3.5) 

For a homogeneous boost, n® =const., this is an exact solution of the equations of 
motion. 

Now, we replace the parameters tq and n® by slowly varying function tqIx) and 
n®(x) of the boundary coordinates x^. Moreover, we promote a, p and the constant 
part of Ai ( which is usually gauged away) to space-time dependent functions. 

ds^ = — v‘^{x)r‘^)dt^ + 2rdtdr — 2r‘^v\x)dt dx® -|- r‘^{dx''Y (3.6) 

( 3 . 7 ) 

A = a(x) (r® — rQ(a;)) dt — a{x)r‘^dr + Ai{x){dx^ — v^{x)dt), (3.8) 

(3.9) 


= fi{x)r 


-6 
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where Ai{x) comes from the constant part of Ai bnt now is replaced by fnnctions 
of x^. The above are no longer a solntion of the eqnations of motion, and we must 
introduce additional correction terms; 


9^iy 9^u T , 

Afj, A^ Oi^ , 

(j) = (j) + if , 

where the background helds A^ and (j) are given by 


-dH) 


(3.10) 

(3.11) 

(3.12) 


4. The first order solution 

In order to obtain the hrst order solution for the hydrodynamic ansatz, we consider 
the derivative expansion. Then, the equations of motion can be treated as ordinary 
differential equations with respect to r, and the correction terms, and (p, can 

be calculated order by order. 

The differential equation can be solved at any given point, and we can take the 
point to be = 0 without loss of generality. The parameters which are replaced by 
slowly varying functions can be expanded around = 0 as 


n*(x) = n*(0) + xf^{d^v^){0) + • • • , 

ro{x) = ro(0) + x^{df,ro){0) H- , 

a(x) = a(0) + x^{d^a){0) + • • • , 
/i(x) = /i(0) + x^{df,fi){0) H-. 


(4.1) 

(4.2) 

(4.3) 

(4.4) 


The derivative expansion of the equations of motion gives the linear differential equa¬ 
tions for the correction terms and ip at the hrst order. The next-to-leading 

terms in ( |4.1|) - (|4~^ ) are at the hrst order and give the source terms in the diherential 
equations. 

Solving the (inhomogeneous) linear diherential equations for the correction terms 
h^u, CLfj. and (p, we obtain the hrst order solution in the derivative expansion. (See 
Appendix ^ for more details.) The integration constants generically modify the 
source terms near the boundary. These contributions can be eliminated by setting 
the integration constants that modify the sources to zero. The hrst order solution 
for the metric is given by 


= —Afdt^ -I- 2rdtdr -|- r‘^{dx^ — v^dt)"^ 

+ ‘^r'^div'^dt^ + AFaij{dx^ — v^dt){dx^ — v^dt) , 
o 


(4.5) 


where cr^ is the shear tensor 


[diV^ + djv") - -dkv’^S, 


i] 1 


iXij — 


(4.6) 





and the function F{r) is given by 


F{r) 



and integration constant is chosen such that F(r) ^ 0 in r —)• oo. 
The hrst order solution of the gauge field is 


( 4 . 7 ) 


A = a(x) 


(r® - ro(a;)) - -r^div\x) 


dt — a{x)r‘^dr + A.i{x){dx'^ — v\x)dt) , (4.8) 


and the dilaton has no correction term, <^9 = 0. Equations of motion implies that 
(p.9|) and its derivatives must be satisfied even after /x and a is replaced by functions. 
The solution above must satisfy the following constraints; 


dta + v'^ditt — adiv\ 

( 4 . 9 ) 

dtTo + v^diTo + ^rQdiv\ 

( 4 . 10 ) 

dtAi + v^djAi -|- AjdiV^ + bar^diro . 

( 4 . 11 ) 


5. Calculation and renormalization of the boundary stress 
tensor 


In this section, we consider the stress-energy tensor on the boundary. In order to 
study the asymptotic behavior, we introduce the vielbein. The leading order term 
of the induced metric is expressed as 

+ r‘^5abe^^el (5.1) 


and 






( 5 . 2 ) 


where — n^Uu, and g^u is given by (|3.6|) and is the normal vector to the 

boundary. On the boundary dr = 0, the vielbein are given by 


T^dx^ = dt , (A^dx^ = dx^ — v°‘dt , (5.3) 

= Vt + v^V^ , = V. . (5.4) 

We also express the gauge field in this frame; 

Aq = v^A^, , Aa = e>^Ai, . (5.5) 

The variation of the action is expressed as 

5S = J + J°5io + >54 + ■ (5-6) 
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Then, the stress-energy tensor and current are given by 

f\ = Sy - , 

jr- = + j“e; . 


(5.7) 

(5.8) 


It should be noticed that the stress-energy tensor is not a symmetric tensor. It 
is related to the ordinary Brown-York tensor as 

= T(by)^ + rAy , (5.9) 

where is the Brown-York tensor, which can be expressed in terms of the extrinsic 

1 


curvature K^j_y as 


rjiflU _ 

(BY) - gvrG 




The current is calculated as 

= 


1 5S 




(5.10) 

( 5 . 11 ) 


V -7 

where is the unit normal to the boundary. 

In order to obtain a regular result, we have to introduce the boundary counter 
terms; 

Sc. = tAt; f d‘xV^ (-7 + . ( 5 . 12 ) 


IGttG 

Then, the stress-energy tensor is obtained as 

1 


o 

o 

1 

1 

00 


T*o = 

1 

1 

CX) 


^0 _ 

i — 

1 

1 

00 

ar 

rjii 

3 ~ 

1 

1 

CX) 

q 

,^5 


3 5 
2 0 


-v^A, + 0(r-^) , 
a 


) + C»(r 


-A + 0{r-*^) , 

^0% - ) + 0{r-'^) 

Zi a 


( 5 . 13 ) 

( 5 . 14 ) 

(5.15) 

(5.16) 


5.1 Energy and momentum conservation 

We now consider the conservation of the stress-energy tensor. Hereafter, we omit the 
factor of r~^ when we discuss on the regular quantities on the boundary. 

In the Newton-Cartan theory, the conservation law takes a slightly different form 
from standard relativistic cases. It cannot be expressed in a unihed form in terms of 
the space-time stress-energy tensor, and we have to introduce the energy vector 
momentum density and stress tensor T^y, which are dehned by 


V = r e 

' /i u' p^a ^ p ’) 


= TPJ& 


‘)(eVe^ 


p 


(5.17) 

( 5 . 18 ) 

( 5 . 19 ) 


- 10 - 



Then, the conservation of energy and momentnm is given by [pi 


38 


- v^{v>^v,). 


(5.20) 

(5.21) 


From the hrst order solntion, ( b.l3|) - (|5.16 ) with C = 1, the energy vector mo¬ 
mentnm vector (1-form) and stress tensor T'^u are given by 


= 


IQttG ° ’ 




T* = 

i6xG . 

(5.22) 

v^ = 

^ A 

IGvrGa * ’ 

(5.23) 

r\ = 

SnG WtiG ° ^ 

(5.24) 


The other components of T^u vanish. We observe that they satisfy the Lifshitz- 
scaling invariance condition 

zS^ = SijTj . 


Using (|5.22|) - (|5.24|) , the energy conservation (|5.20|) becomes 


0 = 


15rP,r„ + 15riv%r, + Sr^d^ - 5rPV„ - 20O.r„)^ - 


(5.25) 


(5.26) 


and its leading order terms give 

5 


(Srpfro + ^rgV^diro + rgdy ) ^ 0, 


(5.27) 


which is the constraint (|4.10|) . The momentnm conservation ( b.21|) becomes 


1 1 ■ 1 1 

0 = ro^jro -dtAi -v^djAi + -AjdiV^ - -dj{rlaij) Ai 

a a a 2 


The leading order terms give 


dt 


dj 


(5.28) 


1 1 1 • 

^r^diVo + -AjdiV^ + -dtAi -F -v^djAi -h A 

a a a 


dt 


dn 


= 0 (5.29) 


which is a combination of ( [4.9|) and (|4.11| ). 

The next-to-leading order terms of the conservation law provide the constraints 
at second order. In order to calculate the solution to second order, we need to 
introduce the derivative expansion of the correction terms, for example. 


dfiu — dfiu + - 


(5.30) 


11 

































where is given by (|3.6|) and is the correction terms which we calculated in 
the previous section. The expansion parameter e is that of the derivative expansion, 
dfy = 0{e). In order to calculate the second order solution, we further introduce 
the second order correction terms hl^u- However, as the correction terms do not 
contribute to the constraint at the hrst order, these second order correction terms 
do not contribute to the constraint at the second order. Therefore, we do not need 
to take hill into account to study the second order constraints. 

The background metric g^y does not consist only of 0{e^) contributions, but also 
contains higher order corrections. The higher order corrections in g^y are included 
in a:-dependent parameters ro, v\ a, and Ai. They can be expanded as 


ro((r) 

= (x) + (a;) H-, 

(5.31) 

v\x) 

= n*^°^(a;) -1- ev^^^^x) + • • • , 

(5.32) 

a{x) 

= a^°^(a;) -|- ea^^^(x) + • • • , 

(5.33) 

Ai{x) 

= (x) + eA\^^ (a:) H-. 

(5.34) 


where etc. are the leading order terms which we studied in the previous section, 
and satishes the constraints at the hrst order. The higher order terms, for exam¬ 
ple do not contribute the hrst order terms, but must satisfy the second order 
constraints. 

After some algebra, the second order constraint equation for the gauge held gives 
0 = . (5.35) 


Together with the hrst order constraint, it can be expressed as 

0 = dta -\- v’^ditt — adiv'’ (5.36) 


This implies that there are no additional terms in this constraint at second order. In 
a similar fashion, from the spatial component of the constraints in Einstein equation 
we obtain 

0 = -divl - -dtAi - -v^djAi - -AjdiV^ + • (5-37) 

a a a 2 

From the temporal component, we obtain 


0 = 15r„8,r„ + 15r,v'd.r„ + - 5r„8 r„ - 20(ftr„) - . (5,38) 


(^.38|) agrees with (|5.26|) , and an appropriate combination of (|5.36|) and (|5.37|) gives 

(Pi. 

From the stress-energy tensor and the current, we can read oh the energy density 
S, charge density n and pressure P as 

3 , 1 _ 1 


8 = 


IGttG 


0 ) 


n = 


IGnGa ’ 


P = 


8nG 


(5.39) 
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In terms of these quantities, the energy flow, momentum density and stress tensor 
are expressed as 


where 


K = 


SttG ° ’ 

In terms of the temperature 


fluid variables and transport coefficients can be expressed as 



(5.40) 

£P - KdiT , 

(5.41) 

nAi , 

(5.42) 


(5.43) 

" = le.G’'" ■ 

(5.44) 

_ 5 2 

“ 47r^° ’ 

(5.45) 


£ = 




and 


IQttG V 5 / 

1 /4;r 


P = 




SttG \ 5 / 


K = 


. —T 
SttG \5 J 


V = 


4:71 ^ \ 

■ —T 

167rG V 5 J 


\ 3/2 


The scaling dimension under the Lifshitz scaling is given as 


(5.46) 


(6.47) 


[£] = 5 , 

= 5 , 

(5.48) 

[n] = 3 , 

IT] = 2 , 

(5.49) 

V] = 1 , 

|.4i] = 1 , 

(5.50) 

[k] = 3 , 

[l] = 3 . 

(5.51) 


In terms of the above fluid variables and transport coefficients, the fluid equations 
take the following form; 


0 = dt£ + v"di£ + {£ + P)diV^ - - diindiT) , (5.52) 

0 = diP + ndtAi + nv^djAi + nAjdiV^ — dj (rjaij) , (5.53) 

0 = dtn + dj{nv^) . (5.54) 

Eq. (|5.53|) can also be expressed as 

diP - dj (rjaij) = (5.55) 

where the current is given by 

= nv^^ . (5.56) 
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V is defined in (^) and the field strength is dehned as = dA with 

A = Ai(dx^ — v^dt) . (5.57) 

To compare, the ordinary non-relativistic flnid eqnations are given by 

0 = dtS + v'-diS + {8 + P)diV^ - - di{KdiT) , (5.58) 

Q = diP + ndtv^ + nv^djV^ — dj {rfaij) , (5.59) 

0 = dtn + dj{nv^) , (5.60) 

where ( |5.58|) gives conservation of energy, the Navier-Stokes eqnation ( b.59|) comes 
from conservation of momentnm, and the continnity eqnation ( |5.6U[ ) implies the con¬ 
servation of mass density. Eqs. (|5.52 ) and ( 5.54 ) agree with the energy conservation 
( b.58|) and continnity eqnation (|5.60|) , respectively. However, ( b.53|) is different from 
the Navier-Stokes eqnation (|5.59|) . In the Newton-Cartan theory, the Milne boost is 
defined snch that the torsionless part of the Newton-Cartan connection is invariant 
nnder the following transformation; 


^ , 


A^A + - -h^^^'V^V^TpdxP 


(5.61) 

(5.62) 


where 14 is the vector which parametrize the Milne boost, is the inverse timelike 
vielbein, is the inverse metric on the time slice, Ap is the gange held in the 

Newton-Cartan theory and P^ = 5'^ — Althongh the inverse timelike vielbein 

wonld be given for the hat spacetime as = 6'^t, h is replaced by = (l,n*) in 
onr conservation eqnation (|5.20|) and (|5.21|) . When we transform to v^, the 

gange held in the Newton-Cartan theory A transforms as 


A ^ A = A + v^dx^ 


1 


v'^dt . 


(5.63) 


If the gange held for vanishes, namely, Ap = 0, we obtain 

A = v^'dx'' — . 

2 

Then, the Navier-Stokes eqnation can be expressed as 

= diP - dj (ijaij) , 


(5.64) 


(5.65) 


where P = dA. Then, ( b.53| ) agrees with the Navier-Stokes eqnation in Newton- 
Cartan theory ( |5.65|) , by identifying the gange held Ap as that in the Newton-Cartan 
theory A. Comparing (|5.57|) with (|5.63|) , we obtain the corresponding Newton- 
Cartan gange held before the boost, A as 


A = (Ai — n*) dx^ — 


v'^Ai — -v^ 1 dt 


(5.66) 


14 






























As is discussed in ||^ , we choose the gauge held A such that 


(5.67) 


where = (1, 0) and P^^' = diag(0,1, • • ■ , 1) is the projection to the spatial direc¬ 
tions. This implies 

1 


A, = . 


z A, 0. 

Then, ( ^.53|) is expressed as 

diP + ndtv'- -f nv^djv'' - dj {Wij) = ■ 


(5.68) 


(5.69) 


This is the Navier-Stokes equation with a non-trivial source term. The interaction 
to external source, comes only from the interaction to the Newton potential 

<h = At since Ai = 0,^’^ and (|5.69| ) can be expressed as 

diP -I- ndtv^ + nv^djV^ — dj (ijaij) = ndi^ . (5.71) 

This expression agrees with the Navier-Stokes equation with the external force T); 

diP + ndtv^ + nv^djv'‘ — dj {r]aij) = Fi , (5.72) 

with gravitational force from the Newton potential T) = ndt^. Thus, ( b.53|) can be 
interpreted as the Navier-Stokes equation in the non-trivial Newton potential; 

1 


$ = A = — 


(5.73) 


The energy conservation ( b.58|) is not affected by the Newton potential since it 
is for internal energy. The conservation of total energy is obtained by appropriate 
combination of (|5.52|) and (|5.53|) , and can be expressed in terms of the Newton 
potential (|5.73|) as 


Q = dt[ S + 2^^^^ — n^] + di 


£ + P + 2^^^^ ~ ~ ~ 


-|- ndt^ . 
(5.74) 

This is consistent with the energy conservation of fluids in the Newton potential; for 
time-independent Newton potential dt^ = 0, the energy conservation is expressed as 


0 = dt { £ + ~ +9i 


£ + P + J n* — rjaijV^ — ndiT 


. (5.75) 


^The Newton-Cartan theory gives Newton gravity for flat Galilean manifold with Fj = 0. The 
Newton gravitational field gi is given by gi = Fu. 

^From (5.67), the Newton potential is expressed in terms of A as 


$ = + -P^'^A^.A, . 


(5.70) 
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(5.76) 


The stress-energy tensor ( ^■13|) -( ^T^) is now interpreted as 


T°o = - -n^ + , 

T*o = — ^ + P — -I- n* -|- rjaijV^ + ndiT , 

= nP , 

T'^j = P5ij — rjaij + nv^v^ . 


(5.77) 

(5.78) 

(5.79) 


This is the ordinary stress-energy tensor for non-relativistic flnids with the Newton 
potential terms. The conservation of total energy (^.74|) and Navier-Stokes eqnation 
( ^.71|) can be expressed in terms of this stress-energy tensor as 


= 0 


(5.80) 


6. The entropy current 


The holographic entropy cnrrent Jg is dehned by the dnal of the {d — l)-dimensional 
volume form on the time slice on the horizon; 


A • • • A dx^’^. 


( 6 . 1 ) 


The entropy current Jg can be expressed in terms of the normal vector as 


\/h 




= 

^ 4G nO ’ 


( 6 . 2 ) 


where the normal vector at the horizon is given by 


rif, = df,S , 


S = r — ro{x) , (6.3) 


to hrst order in the derivative expansion.® 

Then, the holographic entropy current for the Lifshitz space-time with <7 = 4 
and 2 ; = 2 is obtained as 


~ 4G ° 


= 4G’-' 


3„.i 



(6.4) 

(6.5) 


Comparing this expression 
satishes 


with (|5.22|) and (|5.39|) , we hnd that the entropy current 


TJ^ = S>^ + + Pv^ , ( 6 . 6 ) 

®The horizon radius has corrections at higher order in the derivative expansion. Since these 
correction terms appear from the second order, it does not contribute to the entropy current at the 
first order. 
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where the Hawking temperature T is given by 


T = 


47r 


0 ! 


(6.7) 


to hrst order in the derivative expansion. 

We can easily check that the entropy current satishes the second law. The 
divergence of the entropy current is calculated as 

1 






By using (^.37|) and (p.SSP, it becomes 


1 1 

d,.Js = ^Mdirof + -^roaijaij > 0 , 


That is manifestly non-negative. It therefore satishes the second law. 
The entropy density s is given by 

„ _ 7-0 _ J_„3 

® “ 4G " ■ 


The ratio rj/s saturates the KSS bound 


a 

s 


1 

dvr 


( 6 . 8 ) 


(6.9) 


( 6 . 10 ) 


( 6 . 11 ) 


7. General background gauge field 

We have imposed the regularity condition at the horizon r = tq to the gauge held 
and considered the following ansatz; 


A = a{x) (r® — ^ 0 ( 0 ;)) dt — a{x)r^dr + Ai{x){dx'' — v\x)dt) 
This ansatz implies a non-zero Newton potential 

5 = A + if" = -\v ^. 


(7.1) 


(7.2) 


However, as is discussed in Appendix even if we do not impose this condition, 
the singularity appears only at the past horizon. Therefore, we can consider a more 
general background for the boundary theory. 

To implement this, we consider a more general ansatz; 


A = a{x) (r® — rQ(x)) dt — a{x)r‘^dr + At{x)dt -|- Ai{x)dx^ . 

The solution for the correction terms is not modihed by this generalization. 


(7.3) 


A = a(x) 


— 'rQ(x)) — -r^div\x) 


dt — a{x)r‘^dr + At{x)dt -|- Ai{x)dx'^ , (7.4) 
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but the constraint ( |4.11| ) is modified to 


divl = 




(7.6) 


where 


J- = dA. , 


A = Atdt + Aidx^ 


By identifying A to the gauge field in the Newton-Cartan theory as 


A = A = A + v^dx^ — -v'^dt , 


the constraint (|7.5|) becomes 


dirl + -dy + -vwy = . 

a a a 


Taking into account the second order terms we obtain 

diP + ndy + nv^y - dj {rjaij) = . 


(7.6) 


(7.7) 


(7.8) 


(7.9) 


This is in agreement with the non-relativistic Navier-Stokes equation ( 5.59 ) with an 
external gauge field term, which contains the Newton potential in F as = <h. 
The energy conservation (|5.52|) does not change but an appropriate combination with 
(O) gives the conservation of total energy in the presence of the external field; 


0 = + di 

+ ndy + nv^iAj . 


8 + P + — n<h I n* — rjaiy — ndiP 


(7.10) 


Now, there are no constraints on the Newton potential <h and Ai, and hence it is 
arbitrary. 

Next, we consider the stress-energy tensor. It becomes 


T°n — 

1 

1 


SttG P 


1 

1 

0 — 

SttG P 

^0 _ 

i — 

1 

1 

SttG ar 

rjii 

1 

1 

3 

SttG P 


- ^-PA^j + 0{ry , 

+ c){ry , 

A + o{ry , 

- “(A + - ylaij + -h 0{ry 

Qj Zi Qj J 


(7.11) 

(7.12) 

(7.13) 

(7.14) 


Recall that this stress-energy tensor is not defined in the ordinary fashion but includes 
contributions from the current and gauge field. We propose that the above stress- 
energy tensor should be identified with the physical stress-energy tensor as 


rpll _ rplM _l_ JfJ. A _ XII Jp A 

J- 1/ ± 1/ ^ U J^if u i/U -AXp 5 


(7.15) 
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where = nv^. Now, the energy flow and stress tensor deflned by using the 
physical stress-energy tensor are the same to those deflned by (|5.22|) and 
( b.24|) , respectively. However, the momentum density is absorbed into the external 
gauge flelds and becomes 

Vi = 0 . (7.16) 

This is consistent with the Ward identity of a Milne-invariant theory 


= 0 . (7.17) 

The conservation equations in the Newton-Cartan theory with external gauge held 

terms 

, (7.18) 

- V^(TO) + , (7.19) 

agree with the constraints from the bulk equations of motion and become 

^ = dtS + v^diS + {8 + P)diV^ - - di{KdiT) , (7.20) 

0 = diP - dj (rjaij) - , (7.21) 

0 = dtu + dj{nv^) . (7.22) 


Then, these equations agree with the non-relativistic fluid equations if the gauge held 
is given by (Q and A = 0. Since the energy flow and pressure are the same to 
those in the previous sections they satisfy the same thermodynamic relations 


TJ^ = 8^^ + Pv^ 

= 

= + Pv^^ . 


We can further define the following stress-energy tensor, 

fp = + PA, = fpp + 5PPAp . 

By using the identification ([f.TD, the stress-energy tensor is expressed as 

T^o = — { 8 — nd) -I- , 

'' 2 W,i 


pQ = - \^8 + P - + -nv j n* rjaijv^ + udiT , 

T°j = nn* -|- nAi , 

Pj = P5ij — rjaij + nv^v^ + nv^Aj . 


(7.23) 


(7.24) 


(7.25) 

(7.26) 

(7.27) 

(7.28) 








where the Newton potential is identihed as <h = At- The above stress-energy tensor 
contains contributions from the external source helds Afj,. These contributions can 
be subtracted by taking the following combination; 

f• (7.29) 

The physical energy current physical momentum density and physical stress 
tensor are dehned by 




(7.30) 

A = , 


(7.31) 

r^ = T^(e“e:))(eA^) • 


(7.32) 

= diag(0,1,1,1), or equivalently, 

8^ = -T% , Pi = f\ , 

'-f-'i rpi 

' j ~ j ' 

(7.33) 


The above dehnitions of the physical quantities includes the contribution from the 
external held. In order to obtain those without external helds, should be replaced 
by (TT^) ^ 

For ^ = 0, the following conservation equations reproduce the constraints in the 
bulk equations of motion; 


= 0 . 

= 0 , 

and agree with the non-relativistic huid equations; 


Q = dt[8 + -nv"^ ) di 


8 + P + ) n* — rjaijV^ — ndiT 


0 = dt{nv^) -|- dj {PSij + nv^v^ — r^aij) , 
Q = dtn + di{nv") , 


(7.34) 

(7.35) 


(7.36) 

(7.37) 

(7.38) 


After substituting ( [7.38|) to the other two equations, they reproduce with the con¬ 
straints the equations of motion in the bulk, and agree with the non-relativistic huid 
equations. For A ^ 0, the conservation law becomes 


+ J^d^A. = 0 , 


= 0 . 


(7.39) 

(7.40) 


The conservation of total energy and the Navier-Stokes equation are expressed as 


0 = dt { 8 + -nv^ — n<F + di 


8 + P + 2^'^^ ~ 1 


rjaijV^ — K-diT 


+ ndt^ + nv’‘dtAi , 

0 = diP + ndtv^ + nvWjV^ — dj (rjaij) — , 


(7.41) 

(7.42) 


20 










where contains the Newton potential <h = At- 

The energy flow, momentnm density and stress tensor constrncted from the new 
definition of the stress-energy tensor are 

8 ^ = 8 - , 

£:* = 8v^ — nv^v^A/j, — ndiT , 

Vi = nAi , 

V j PSij ^o'ij . 

Then, the thermodynamic relation is expressed as 

TJ^ = 8>^ + jlJ^^ + Pv>^ 

= A + Pv>^ 

= + Pv^ , ( 7 . 47 ) 

where the chemical potential fi is given by 

fi = • ( 7 - 48 ) 


(7.43) 

(7.44) 

( 7 . 45 ) 

( 7 . 46 ) 


The flnid in the Newton-Cartan theory was stndied in [^. The relation between 
the variables in this paper and those in ||^ is presented in Appendix The flnid 
obtained in this paper is the same to that in |^, but some transport coefficients are 
absent here. There are 6 transport coefficients in and 3 of them are related to 
parity-odd terms. Such terms are not present here as we ignored all potential CP- 
odd terms in the bulk effective action. Adding the gravitational, gauge and mixed 
Chern-Simons terms we expect that such CP-odd terms will be generated but we 
leave this for future work. 

The bulk viscosity vanishes here because of the Lifshitz scaling symmetry. Fi¬ 
nally the remaining two transport coefficients, the shear viscosity and heat conduc¬ 
tivity appear in our model and take concrete values associated to the bulk theory in 
question. Furthermore, we study only the flat background for the geometry on the 
boundary, and e^, and hence, related quantities such as torsion, do not appear in 
this paper. 


8. The case of general z 

So far, we have focused on the case of 2 ; = 2. In this section, we consider the 
hydrodynamic ansatz for z > 1. 
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As in the case of z = 2, we first replace the parameters by the slowly varying 
functions. For general we start from the following ansatz; 


ds^ = — v^{x)r‘^)dt^ + 2r^ ^dtdr — 2r'^v\x)dt dx^ + r‘^{dx''Y 

/ = 1 


rZ+3 


A = [a{x) — rQ'^^{x)) — Ai{x)v\x)] dt — a{x)r‘^dr + Ai{x)dx'', 

= n{x)r~^, 


( 8 . 1 ) 

( 8 . 2 ) 

(8.3) 

(8.4) 


and calculate the correction terms by using the derivative expansion. 

As for z = 2, we introduce the correction terms and solve the linear differential 
equations. The integration constant can be hxed by the requirement for the asymp¬ 
totic behavior and the regularity at the horizon. The hrst order solution becomes 

ds^ = -r^^fde + 2A-^dtdr + r\dA - v^dtf 

+ ‘^AdiV^dt^ + 2F^{r)dirQdt{dx^ — v^dt) — r‘^Fi{r)aij{dx^ — v^dt){dx^ — v^dt) 

(8.5) 


and 


where 


A = a I 


diV^ ) dt — ar'^dr + F 2 {r)diro{dx^ — v^dt) , (8.6) 


FAr) = / dr- 


T — Tn 


+3 _ 


lyt tiyt A, “p Kj 


(8.7) 


F 2 {f) = (2{z — l)r^’''^ ~ ~ 5)rQ’''^) 

{z -f 3)ar^rg’''^ [10(^ “ l)r^+^rQ -|- z{z -f 3)r®rg — (z — 5)(z — 2)r^+^] 


X / dr- 


2{z - l)(r^+3 - ro+^)[2(x - l)r^+3 - {z - 5)ro+^]2 


r 7 ^ 2(x- 1) 

Fsir) = - I ^- F2{r) , 


a 


( 8 . 8 ) 

(8.9) 


where the integration constants are determined by the asymptotic behavior at the 
boundary, r —)■ cxd, 


Fi(r) = 0{r-^) , F2{r) = , Fs{r) = , 


( 8 . 10 ) 


for 1 < X < 2. For x > 2, they have different asymptotic behavior but the integration 
constants can be taken to be the analytic continuation of those for 1 < x < 2. 
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The functions ro, u*, a and Ai in the solution must satisfy the following con¬ 
straints; 

: dta + v'^ditt — adiV 
1 


0 


These constraints agree with the leading order fluid eqi 
In order to obtain the regular results, we introduce 


( 8 . 11 ) 

( 8 . 12 ) 

(8.13) 


0 = dtro + v'diTo + ^rodiv\ 

0 = dtAi + v^djAi + AjdiV^ + • 

2 ( 2 : - 1) 

leading order fluid equations. 

f results, we introduce the following counter terms; 

. (8.14) 




Then, the stress-energy 
1 1 


^0 _ _ 

° 87rG' 

fr _ L 

° 87rG' r^+3 


-(5 + ^)+ 2 

tensor is calculated as 

■vA, 






^0 _ 

i — 


\ 

C)(^-h+4)) 

■A: 


a 

+ A±lLrfa^r„ - —+ A,v 

Qj Zi 


(8.15) 


. - nr 

- 1 ) ' 




87rG ar^+3 


_^2+3x. 

O ' 0 


0 *^0 2 A 


^ - .i 


rj^t ___ 

^ 87rG'r^+3 1^2 “ " 'i - ^ a 

the above expression of the stress-energy 

z-1 


From 
ables as 


f — ^z+3 „ _ ^ p _ ^z+'d 

IGttG ° ’ IGirGa ’ IGttG ° 

the energy density and pressure satishes the scaling invariant condition 

zS = {d- 1)P. 

momentum density and stress tensor are 

T° = T , 

= SA - KdiT , 

Vi = nAi , 

P j Pdij ^(Zij , 


(8.16) 
(8.17) 

■'ti'U, j + 0(r-<'+‘‘>) . (8.18) 

tensor, we identify the fluid vari- 

(8.19) 


p -2 + 3 


as 


where 


Here, 


The energy flow. 


I expressed 


S{z-l)G 


'^0 




( 8 . 20 ) 

( 8 . 21 ) 

( 8 . 22 ) 

(8.23) 

(8.24) 

(8.25) 
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The fluid variables and transport coefficients can be expressed in terms of the tem¬ 
perature 


T 


z + 3 
An 


0 > 


(8.26) 


as 


£ = 


An 


IQnG \ 2 : -f 3 


■T 


£+3 


p = 


IQnG V -f 3 


An 


■T 


2 + 3 


and 


K = 


8(2- 1)G V^ + 3 


An 


-T 


z + l 


h = 


An 


IQnG \ 2 -|- 3 


-T 


(8.27) 


(8.28) 


The scaling dimension under the Lifshitz scaling would be 


[^] 

— z d — 1 , 

+ 

1 

(8.29) 

N 

= d — 1 , 

[T]=z, 

(8.30) 

V] 

= z-l , 

= 1 , 

(8.31) 

[/.] 

= z d — 3 , 

[ri]=d-l . 

(8.32) 


The conservation law can be understood as the following fluid equations, 

Q = dt£ + Pdi£ + {£ + P)div'- - ^Wij(^ij - di {ndiT) , 

0 = diP + ndtAi + nv^djAi + nAjdiV^ — dj (rjaij) , 

0 = dfU -I- dj{nv^) . 

For general z, the entropy current is calculated as 

1 , 


- 4g''5 • 

P — J_.r 3 i_ ^ 

~ AG 8(2 - 1)G 


r^diro . 


They satisfy the the thermodynamic relation; 


TJi^ = £f^ + Pv^ = + Pv^ ■ 


In terms of the temperature, entropy density is expressed as 


s-A- — (-^T 

Ag[z + 3-^ 


(8.33) 

(8.34) 

(8.35) 


(8.36) 

(8.37) 


(8.38) 


(8.39) 
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9. Results, interpretation and outlook 

In this paper, we have investigated the fluid/gravity correspondence for Lifshitz 
invariant theories. We have considered the Einstein-Maxwell-dilaton theory which 
has the Lifshitz geometry as a solution. The gauge held and dilaton break however 
this scaling symmetry mildly and this breaking is characterized by a non-trivial 
conduction exponent (. The geometry which describes hydrodynamics is constructed 
from the black hole geometry. We have used the Eddington-Finkelstein coordinates 
to impose the regularity condition at the future horizon, and we have then introduced 
a Galilean boost. We have replaced the boost parameter and horizon radius by slowly 
varying functions of the space and time coordinates. We have done the same for the 
parameters appearing in the gauge held and the dilaton. The solution is derived to 
hrst order but the constraints are calculated to the second order. 

The relation between the bulk constraint equations and the conservation laws on 
the boundary implies that the boundary theory satishes the conservation laws in the 
Newton-Cartan theory. The conservation laws have the form of huid equations but 
in a non-trivial Newton potential, if we impose the ordinary regularity condition of 
the gauge held. At = 0 at the horizon. However, the geometry which describes the 
huid is generically not regular at the past horizon, even if the regularity condition 
at the past horizon is imposed on the gauge held. Imposing the regularity condition 
on the gauge held only at the future horizon, we obtain the huid equations with 
arbitrary Newton-Cartan gauge held, which contains an arbitrary Newton potential. 

The huid has the following properties; 

• The stress-energy tensor has a form similar to that of non-relativistic huids, 
and is expressed in terms of the huid variables: the velocity held n*, the energy 
density S, the pressure P and the charge density n. It also contains the exter¬ 
nal gauge held A. At hrst order, it has as transport coefficients the thermal 
conductivity k and shear viscosity 77 while the bulk viscosity is zero due to the 
Lifshitz scaling symmetry. All the above variables and constants are functions 
of temperature T. 

• The (particle number) density appears associated to the external gauge held. 
The stress-energy tensor agrees with that of non-relativistic huids except for 
the terms where the density appears. The terms which contain the density are 
diherent from those of a non-relativistic huid. In particular, the momentum 
density is proportional to Ai and vanishes for Ai = 0. 

• The conservation law of the stress-energy tensor is not given in the form of 
its covariant derivative, but is that appropriate to the Newton-Cartan theory. 
It is expressed in terms of the energy how momentum density Vt and 
stress tensor T*j. The stress-energy tensor satishes the Ward identity of the 
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scale invariance zS = {d— 1)P without any modification. In terms of the fluid 
variables, the conservation law takes the form of the fluid equations: the energy 
conservation, continuity equation and the Navier-Stokes equation. 

• The fluid equations are similar to the non-relativistic fluid equations. The 
continuity equation and energy conservation equation agree with that for the 
non-relativistic fluids. However, the Navier-Stokes equation is different from 
that for ordinary non-relativistic fluids. It does not contain the velocity field 
except the terms appearing in the shear tensor. It contains the coupling to the 
external gauge field instead. 

• The absence of the velocity fields in our Navier-Stokes equation can be ex¬ 
plained as follows. For an ordinary fluid, the pressure is much smaller than 
the energy density. In our case, the pressure is of the same order as the en¬ 
ergy density because of the Lifshitz scaling symmetry. If it is interpreted as 
a non-relativistic limit of a relativistic fluid with Lifshitz scaling symmetry, 
the contribution from pressure is much larger than the one coming from the 
velocity fields. This is the reason that some terms with velocity fields in the 
Navier-Stokes equation are absent. These terms are replaced by the external 
source which is identified with the gauge field in the Newton-Cartan theory. It 
should be noted that these terms does not vanish if the non-relativistic limit is 
imposed by the light-cone dimensional reduction. 

• We may do some redefinitions and by identifying the gauge field A to the 
gauge field in Newton-Cartan theory, which takes the form of ^ = A + v^dx^ — 
^v'^dt after introducing the boost, our Navier-Stokes equation agrees with the 
standard non-relativistic Navier-Stokes equation. 

• Since the gauge field in the Newton-Cartan theory is a generalization of the 
Newtonian gravity theory, a general external gauge field gives a fluid in a non¬ 
trivial gravitational potential. The Newton potential appears in the gauge field 
as <h = At- The Navier-Stokes equation we obtain agrees with the ordinary 
Navier-Stokes equation in the presence of an external (Ngravitational) force. 
The gauge field does not contribute to the conservation of (internal) energy 
density S. The conservation of total energy can be obtained from the conser¬ 
vation of S and the Navier-Stokes equation and agrees with that for ordinary 
non-relativistic fluid in a non-trivial Newton potential. 

• The entropy density is defined in terms of the horizon area and satisfies the local 
thermodynamic relation with energy density and pressure. The divergence of 
the entropy current is non-negative, which is consistent with the second law. 
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The form of the fluid equations is independent of the Lifshitz exponent 2 ; as well 
as of the conduction exponent (. The scaling exponent dependence appears 
only in the constitutive relations. 


There are several obvious interesting questions that remain unanswered by our 
work. The hrst is the extension of our results to hydrodynamics in the presence of 
hyperscaling violation in the metric. This is under way. A naive guess would be that 
the hydrodynamics would be a dimensional reduction of the one found here, along the 
lines described in i0.a- In particular in it was shown that Lifshitz solutions 
with hyperscaling violation can be obtained as suitable dimensional reductions of 
higher dimensional Lifshitz invariant theories without hyperscaling violation. The 
associated reduction of the hydrodynamics will provide equations similar to the ones 
here but with a non-zero bulk viscosity. This needs to be verihed. 

A further extension involves Lifshitz geometries with broken U(l) symmetry. 


This is actively pursued in |55 


An interesting question in relation to the above is: what is the appropriate 
hydrodynamics for QFTs that are RG Flows that interpolate between relativistic and 
non-relativistic theories. To motivate the answer to this question, we consider hrst 
non-Lorentz invariant (but rotationally invariant) hows between Lorentz invariant 
hxed points^, 0. but where the velocity of light in the IR is diherent for that in 
the UV. In such a case, the hydrodynamics of this theory, is relativistic, but with a 
speed of light that is temperature dependent. 

This example suggests that in an (Lorentz-violating) RG how from a GFT (with 
a un unbroken U(l) symmetry that is used to drive the breaking of Lorentz invari¬ 
ance) to an IR non-relativistic scaling (rotational invariant) geometry at an arbitrary 
temperature, the hydrodynamics will be again of the relativistic form (but with a 
general equation of state) and with a speed of light c(T) that is again temperature 
dependent. In the IR, c(T 0) = 00 and the hydrodynamics reduces to the one 
found here with the U(l) symmetry becoming the mass-related symmetry. This is 
nothing else than the standard non-relativistic limit® of the the relativistic hydro¬ 
dynamics while all thermodynamic functions and transport coefficients are smooth 
functions of T (If no phase transition exists at hnite T). Otherwise they follow the 
standard behavior at phase transitions. 

A more general breaking of Lorentz invariance during a RG flow must involve 
higher form helds of tensors in the bulk, or a multitude of vector helds and the 
details of the RG flow become complicated. It is important that such flows are 
analyzed as they hold the key to understanding general non-relativistic flows as well 
as generalized hydrodynamics of the associated theories. 


’^The fact that the speed off light can vary on branes was pointed out first in Q. 

®This is expected to happen along the lines presented in although this needs to be verified. 






Finally a more detailed study of the above issues in the absence of U(l) symmetry 
is necessary. 
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APPENDIX 


A. Notations 


Since in the topic treated in this paper there are a lots of variables involved and a 
lot of redehnition we present here, a list of the variables and their dehnitions with 
comments when necessary. We are present a translation dictionary with the variables 
used by Jensen, and Hartong et ah, [p8[ . 


Variables defined on the gravity side 


fh Boost parameter introduced into the (static) black hole geometry in 


Tq: the horizon radius which is dehned in (|2.12|) . 


a: the coefficient of the ^ term of the gauge held, introduced in (p.6|) . 
Note that a ^ a^dx^. 


Afj_: the constant part of the gauge held, which is dehned in (|3.8|) . This corre¬ 


sponds to the Milne boost-invariant gauge held in the Newton-Cartan theory 
Afj_. In this paper, it is sometimes expressed as the 1-form A = A^dx^. 

v^\ the timelike inverse vielbein on the boundary which is dehned up to the 
factor or equivalently, dehned in (|5.2|) and (|5.4|) and given by = 

(l,r;*) in our model. This corresponds to the Milne boost-invariant timelike 
inverse vielbein, or equivalently, the velocity vector held of the huid. 


• r^: the timelike vielbein on the boundary which is dehned up to the factor of 
r^^f, or equivalently, dehned in ( b.lD and (|5.3|) and given by = (1, 0) in our 
model. This corresponds to the timelike unit normal which dehnes the time 
direction in the Newton-Cartan theory. It is automatically invariant under the 
Milne boost. 


• e“: the spacelike vielbein on the boundary which is dehned up to the factor of 

r^, or equivalently, dehned in (^.1|) and (|5.3|) , and given by e'^dx^ = — v°'dt 

in our model. This corresponds to the Milne boost-invariant spacelike vielbein. 

• e^\ the spacelike inverse vielbein on the boundary which is dehned up to the 
factor of r“^, or equivalently, dehned in ( b.2|) and (|5.4|) , and given by = 
da in our model. This corresponds to the spacelike inverse vielbein. It is 
automatically invariant under the Milne boost. 


Variables in the (boundary) Newton-Cartan theory 
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hpv 

Ap 


Ad 

Ap 

n 

$ 


Jensen 32 

np 



A, 


- 

A, 

P 



HKO 38 


—v^ 


-Mp 

—v^ 

Ad 

(-4,0) 


-Mt 



Table 1: Relation between the variables in this paper and those in |^2| and |^. Cor¬ 
respondence to [^] is read off from the relation of the metric and gauge field in gravity 
side. In |^, is defined such that spacial components of the associated Milne invariant 
gauge field vanish. We can also define such velocity field in our notation, If we 

identify this combination to —v^ in [^ ], we obtain another correspondence. In this case, 
no variables in |^] correspond to the Milne boost invariants in this paper. 


v^: the timelike inverse vielbein defined by = (1,0). This notation is in¬ 


troduced in (|5.61|) to discuss the Milne boost. The timelike inverse vielbein is 
not invariant under the Milne boost but it is covariant. It is sometimes called 
“velocity” but must be distinguished from the velocity of the fluid. 

the induced metric on the time-slice. It is given by = diag(0,1,1,1) 
for a 4-dim space-time. It is not invariant under the Milne boost but it is 
covariant. We dehne this here in order to express in covariant form. 


e“: the spacelike vielbein, which is introduced in (|7.31|) - (|7.32|) . It is given by 


e“ = diag(0,1,1,1) for a 4-dim space-time. It satishes = e“e“. 


e^: the spacelike vielbein, which is introduced in (|7.31|) - ([7.32|) . It is given by 
= diag(0,1,1,1) for a 4-dim space-time. Since it is invariant under a Milne 
boost, it is equal to e^. 

An'- the gauge held in the Newton-Cartan theory. This notation is introduced 


in ( b.62|) to discuss the Milne boost. It is not invariant under the Milne boost 
but it is CO variant. 


A^: the Milne boost invariant combination for the gauge held. It is dehned 


m 


(|5.63|) , A = A + v^dx^ — or equivalently, A = A + — 


^v^v'''h^nTpdxP. In the torsion-free case, (h^, Ap) and (h^, Ap) give the same 
Christohel symbols. They are related to each other by the Milne boost. 


VEVs and fluid variables 


• the boundary stress-energy tensor without counter terms. This is dehned 

in (|5.7|) by putting the boundary at hnite r. 
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the boundary stress-energy tensor with counter terms, or its regular part 
in the section in which we are discussing only on the boundary fluid. This is 


introduced in ( |5.13| )- (|5.16|) 


It is deflned by 

used to deflne Milne boost-invariants for a general background of A. 


d^yJ^Ap. This definition is 


T^y\ deflned by ( [f.24| ), T^y = T>^y + S^yJ^Ap. T^y — J^Ay gives the stan¬ 
dard non-relativistic fluids’ stress-energy tensor. It consists of the physical 
energy vector, physical momentum density, and physical stress tensor. They 
are different from the Milne boost-invariants. 


• J^\ the current without counter terms deflned in (^.8|) . It is regular even 
without the counter terms. It corresponds to the mass current. 

• deflned by —T^yv'^ in ( b.l7|) . It corresponds to the (Milne boost invariant) 
energy vector. In section 7, T^y is replaced by T^^y. 

• Vp\ deflned by T^yTp&^^e^^ in ( |5.18|) . It corresponds to the (Milne boost- 
invariant) momentum density. In section 7, T^y is replaced by T^y. It is 
different from the physical momentum density which is not invariant under 
the Milne boost. 


• T^y'. deflned by T'^(e“e(()(e^e()) in ( b.l9|) . It corresponds to the (Milne boost 
invariant) stress tensor. In section 7, T^y is replaced by T^y. 

• deflned by —T^^yv'" in ( [7.3(J| ). It corresponds to the physical energy vector, 
which contains a contribution from the mass density. 


Vp\ deflned by T^yTpe'^e‘1^ in ([7.31|) . It corresponds to the physical momentum 
density, which contains a contribution from the mass density. 


T^y: deflned by T^(e“e(()(eje^) in (|7.32|) . It corresponds to the physical stress 
tensor, which contains a contribution from the mass density. 


• S: (Milne boost invariant) energy density, or equivalently, internal energy den¬ 
sity. It is introduced in ( p.39|) and equals £^. 

• P: The pressure. It can be read off from the stress tensor and introduced in 

(Hi. 

• n\ it is deflned by 1/a. It is introduced in ( |5.39[ ). It corresponds to the particle 
number density, or equivalently, the mass density. 


• T: the temperature. It can be calculated as the Hawking temperature of the 
black hole ( |2.13|) . 
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Table 2: Correspondence of the fluid variables in this paper and those in |32]. In |32], the 
fluid variables do not contain contributions from the external source. To be exact, in 
the definition of and should be replaced by — J^Au, since they contain 

contributions from the external source. 


Jg-. the entropy current, which is dehned from the volume form in ( p.2|) on the 
time-slice at the horizon. 


s: the entropy density, which is dehned in (|6.10|) . It equals to in the non- 
relativistic case. 


r]-. the shear viscosity. It is introduced in ( ^.43|) and can be read off from the 
stress-energy tensor, or the huid equations. 


k: the heat conductivity. It is introduced in (|5.41|) , and can be read off from 
the stress-energy tensor, or huid equations. 


$: The Newton potential, which is <h = At- It is introduced around (5.71). 


Basic fields and constants in gravity side 


the metric. It is introduced in (^T|) 


A: the gauge held. It is introduced in (|2.1|) . 


(j)-. the dilaton. It is introduced in (p.ll). 

G: The Newton constant. It is introduced in ( p.ll) . 

A: the gauge held with local lorentz indices. It is dehned in (|5.5| ). 
J: the current with local lorentz indices. It is dehned in (|5.6|) . 

the correction terms for the metric. It is introduced in ( p.lOl) . 


a^: the correction terms for the gauge held. It is introduced in ( |3.11|) 
(p: the correction terms for the dilaton. It is introduced in ( p.l2| ). 
R^^-. The Ricci tensor. It appears hrst in (|2.1|) . 


A: the cosmological constant. It appears hrst in (2.1). 
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the energy-momentum tensor in the bulk. It appears in (p.ll|) . 


Tl^^yy The Brown-York tensor which is defined by (^.101) , 


n^\ the normal vector to the boundary or horizon. It appears first in ( b.llj) . 


the extrinsic curvature on the boundary. It appears in (^.10|) . To be 
precise, it is defined on constant but finite r surface. 


the induced metric on the boundary. It is introduced in ( b-ll) . To be 
precise, it is defined on constant but finite r surface. 


B. Calculation of the equations of motion at first order z = 2 

In this section, we calculate the first order solution in derivative expansion around 
(^^ - (|3.9|) . The coordinate can be chosen such that v\x) = 0 at any given point, 
therefore we may take n*(0) = 0 without loss of generality. Although we will work 
in n*(0) = 0 coordinates, the solution for n*(0) ^ 0 can be obtained by boosting the 
solution uniformly. From now on we work at the point = 0, but we omit the 
subscript “(0)”, hereafter. Here, we take the following gauge conditions; 

5'rr = 0 , oc , Tr[^"^/i] = 0 , ttr = 0 , (B.l) 

where = (1,n*). 

The correction terms can be classified by using the 50(3) symmetry along the 
spatial directions. The equations of motion are separated into that for scalar (sound 
mode) , vector mode and tensor mode. The sound mode consists of the following 
components 


htt h\ (Xi ip , (B.2) 

and vector mode 


hti 


Qji 


1 


(B.3) 


and the tensor mode is the traceless part of the metric 

K, . (B.4) 

To simplify the differential equations, we redefine the correction terms for the 
metric 


htt 


htr 


hti 


hij 
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as follows 


9tt = r^i-f + hit), 

h 

9tr 

gti = r‘^hti 

9ij ^ h^xx^ij “ 1 “ ^ij) 


(B.5) 

(B.6) 

(B,7) 

(B.8) 


where hxx is the trace part and hij is the traceless part in x*-components. The gange 
condition gives additional constraints; 


‘^htr T ^hxx 0 


(B.9) 


and the (r, r)- and (r, i)-components of the correction terms must vanish. We define 


hi = -hx 


(B.IO) 


The other correction terms and ip are similar to the definitions ( |3.11|) and (|3.12|) , 
but Qr is eliminated by the gauge condition. 

B.l The sound mode 

Some components of the equations of motion do not become differential equations for 
the correction terms, but give the constraints on the parameters. From the Einstein 
equation we obtain 

where is the normal vector and is the induced metric on r =const. surfaces. 
In fact the above equation contains no correction terms. For the sound mode, by 
contracting the ( [B.ll| ) with = (l,n*), we obtain the following constraint; 

4 

0 = 7 ^ {Sarodiv’’ + 15adtro + 2rodtp) + {adiV^ - dtp) r , (B.12) 


where 


z + d — 1 


- 1/2 


■ (b.13) 

Equation ( |B.12|) must be satished for arbitrary r, and hence the hrst and second 
terms must vanish independently. 


Sarodiv" + IbadtVo + 2rodtp = 0 , 
adiV^ — dtp = 0 . 


(B.14) 

(B.15) 


From ( p.3|) , the r-component does not contain correction terms and gives another 
constraint; 

—adiv'' + 2dtp — dta = 0 . (B.16) 
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After substituting the above constraints to the equations of motions for the sound 
mode, we obtain the following differential equations; 

0 = ) _ Q^/Qr^}l[{r) + 6V6rQh[{r) + V6r^h[^{r) + h\/Qr^htt{r) 

Qj 

+ r®(^"(r) — rrQ(p"{r) + 6r^(p'{r) — rQip\r) + 30r"^93(r) — 2\/6r^c?jU® (B-17) 


0 = ra'l{r) — 4aj(r) + 30ar®h'^(r) + 5\/6aor®<y9'(r) (B.18) 

0 = —12ar‘^div'' + 8a^(r) — 72ar^h[{r) — 18arQh[{r) — 240ar^hi(r) 

+ ?>ar^h'l^{r) + 30ar^h^^(r) + 60ar^hti(r) + 20\/6ar"^9?(r) (B.19) 

0 = 3a (—4r^c?jU* — 6 (4r® + r^) h\(r) + r^h”^{r) + lQr^h[^{r)) 

+ 8a[{r) + 20\/6ar^(^(r) — 240ar^hi(r) + 60ar^htt(r) (B.20) 

0 = -12h; + V6ip' - Srh'l (B.21) 


ar^ r 

+ 8r^h'^^{r) + lbr‘^hu{r) — 5\/6r^<y9(r) — 69*^* (B.22) 


where the first equation originates from (|2.4|) and the second from the f-component 
of (|2.3|) . The third to hfth equations are the (t,r)- and (r, r)-components of 

(p.2|) . The last equation is the trace part of the spatial component of (|2.2|) . The above 
equations are not independent but an appropriate combination gives the constraints, 
which we have already imposed, and hence becomes trivial. 

We first impose the constraints to the parameters and then solve the differential 
equations. The solution for the sound modes is 


htt — 


2^/6 


1 - 


+ TT^diV^ - 


+ 1 


^lfe(r) + 


3r2 

1 


4+41 hi' 


( 0 ) 


^ ) h'.{r) 


hi = + ft2(r) 

+ \/6r®0® — 15r®h2(r' 


at = ai 


= y«) _ +3 


-ar^diV^ 

3 


h2(r) + \ ^rh'.^{r) 


where and af'^ are integration constants and ^2 is given by 


4°) 


hiir) = * 7*'y\ Fi 
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5 1 


54 


+ 


3-7 5 

5.0.^ 




+ 444'’.' ( 4- 


5^4 


3,3 


2 2 ’ 2 ' 2 
;, A(i 3_^), T(i3 + yT^)\ 


)’ 


(B.23) 

(B.24) 

(B.25) 

(B.26) 


(B.27) 
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where pFq and are hypergeometric function and Meijer G-function, respectively 
and Cl and C 2 are integration constants. 


B.2 Vector mode 

As for the sound modes, spatial component of ( |B.lip gives a constraint 

0 = dtAi + AjdiV^ + bar^diVo + (r^ - rl)di{p - a) . (B.28) 

Since this constraint must be satished at arbitrary r, the hrst and second terms must 
vanish independently, 

dtAi + AjdiV^ + bar^diVo = 0 , (B.29) 

di{p — a) = 0 . (B.30) 

Then, the equations of motion for the vector modes are 

0 = r (—5rrg {badiVo + rodta) + (r® — Tq) a-{r) + baAh[t{r)) + (Trg — 2r®) a-(r) 

(B.31) 

0 = (—2 (Sargdirg + djO (rg — r^)) + ar^ (rg — r®) h"-(r) + 4ar (rg — r^) h[t{r)) 

- 2 (r= - a'(r) (B.32) 

0 = 20- + Aar^h'^t + aAh"i — 2Cdip, (B.33) 

where the hrst equation is the x'-component of (|2.3|) and the others are the (t,x*)- 
and (r, x*)-components of (|2.2|) , respectively. 


The solution for the vector modes is 


h.,= ft 


lyt'-i 


2 

ti + -l^ditt 

a 3a 




dr ai(r). 


(B.34) 

(B.35) 


where and a'P are integration constants. The function ai is given by 
= “3(r5 -r5)(2r5 + 3r5)2 
and is expanded around r = rg as 


ar‘^[r^ — rg) — dr'^Gg + r^(8r^° + 18rQ°)djaj 

(B.36) 


3G3 - 2 brldia 

“ 375r7(r - rg) ^ ' 

In order for the solution to be regular at r = rg, we we must take 

26 


C?. = yr^diO 


and for the rest we obtain 


hti = 


( 0 ) ^5 


a- r 


a r-^ 
M , .(0) 


tti = ( r® + |ro ) + ^r'^ditt 


(B.37) 

(B.38) 

(B.39) 

(B.40) 
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B.3 Tensor mode 

The equation of motion for the tensor mode is given by 


0 = 2 (— 6 r® + rQ)h[j + 2 (—r® + rrl)h'lj — + diV^) 

There are no constraints for the tensor mode. 

The solution is 


hjn CT- 


r'^dr 




+ ^4 


dr 


r(r5 - r^) 


where 






= {diV^ + djv") - -dkV^d, 




Regularity at tq implies C 4 = Tq. We dually obtain 

f {r^ — rg)(ir 


C. Calculation of the stress-energy tensor 


(B.41) 


(B.42) 

(B.43) 


(B.44) 


Here, we calculate the stress-energy tensor. For the solution ([4.5|) , the Brown-York 
tensor is obtained as 


87rGT(BY)% = 3 - , 

87rGT(BY)*o = ^ (-4^0^* + ^^r|) + C(r-^) , 

87rGT(BY)°* = + 0{r-^^) , 

87rGT(BY)'i = 4% + > 

to the hrst order in the derivative expansion. 

The current is given by 


16vrGJ° = ^ + + 0{r , 

ar^ 6ar‘ 

IQTiGf = + C(r-^') . 

ar^ ar^o ^ ^ ’ 


From ( b.9D we can obtain the stress-energy tensor as 


(C.l) 

(C.2) 

(C.3) 

(C.4) 


(C.5) 

(C. 6 ) 


87rGf % = 4+4 (-2r^ - -uU* 

r^ \ a 

SirGf^ = + C>(r-6) , 

nro 


^ + C(r- 6 ) 
+ C(r- 6 ) , 


StiGT^j 


4:Sij + 


1 



-Tna, 


o'-’v 


+ -v^Aj + 0 (r-®) 


(C.7) 

(C. 8 ) 

(C.9) 

(C.IO) 
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Since the volume form on the boundary behaves as 




(C.ll) 


the terms of become hnite. Those of 0{r^) diverge at the boundary, r —)■ oo, 

but can be subtracted by introducing the boundary cosmological constant term. We 
can further introduce a boundary counter term proportional to 


Sr,= 


IGttG 


(^-8 + C + , ( C . 12 ) 


This induces a counter term for the stress-energy tensor 

8nGT,Au = -(4 - C)SA + + 0{r-^) . 

Subtracting the counterterm, the renormalized stress-energy tensor becomes 


(C.13) 




(C.14) 


r‘o = + 0(r-») , (C.15) 


From the equations above we can read the energy density S and pressure P 


IGttG ° ’ 


D ^ 4" ^ 5 


(C.18) 


The coefficient G in the counter term can be hxed by the regularity condition 
for the operator dual to the dilaton scalar 0. The vacuum expectation value of the 
operator is given by 

= = (C.19) 

a /—7 00 


For our first order solution, this expectation value is calculated as 


o^ = 


IGtiG \ V 2 r® 




(C.20) 


The regular term is at 0{r while the counter terms becomes 


at = 

167rG2 ^ ^ 167rG 


(^VgG Py/GG^^-^ +0{r-^) . (C.21) 
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The renormalized expectation value is 


o.= 


IQttG 


\/6 (l-C') 


5(l-2C)^|+0(r-») 


(C.22) 


To obtain a finite we must take C = 1, as the first term is divergent at the 
boundary. We obtain 


7^0 _ 

1 

1 

0 

SttG P 

n = 

1 

1 

SttG a 

rriO _ 

1 

1 

i 

SttG ar 

rni 

3 ~ 

1 

1 

SttG a 


+ 0 (r-^) , 

+ ]^dirl - ) + 0{r~^) 


Ai + 0(r-^) , 

^o^ij - \rlaij + ^v^Aj ) + 0{r-^) , 


and 




3r5 


+ C>(r 


- 6 ^ 


IGttG V 2 

Moreover the renormalized stress tensor (|5.13|)-(|5.1(j|) satisfies 


(C.23) 

(C.24) 

(C.25) 

(C.26) 

(C.27) 


z£ = {d-l)P . 


(C.28) 


D. More on counter terms 


We can also consider higher order terms of for the counter terms. Due to the 
constraint on a and p, always appears with the factor of e^^/2. Then, general 
counter terms is expressed as 


S 


Ct 



Q 


rv ^ / 


(D.l) 


where the boundary cosmological constant term is fixed such that the boundary 
stress-energy tensor becomes finite. Then, the stress-energy tensor becomes 


o 

o 

1 

1 

SttG P 

T\ = 

1 

1 

SttG P 

rn{) _ 

i 

1 

1 

SttG ar' 

rni 

1 

1 

^ j ~ 

SttG r® 


{-"A) 


- -vU, 

a 


+ 0{r-^) , 




-A + 0(r-^) , 


1 + - ^r^aij + AAj 


+ 0(r 


- 6 \ 


(D.2) 

(D.3) 

(D.4) 

(D.5) 
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The dual operator to the dilaton becomes 


O0 


1 

IbvrG 





+ 0(r-®) . 


In order to regularize this, the coefficient must satisfy 


Co — 1 . 
a 

Then, the stress-energy tensor is same to (|5.13|) - (|5.16|) . 


(D. 6 ) 


(D.7) 


E. Regularity conditions of the gauge field at the horizon 

If the guage held has non-zero At at the horizon, it becomes singular at the 
horizon. It can be easily seen by taking the wick rotation and by considering the 
Polyakov loop wrapping on the time circle. Although the horizon is a point in the 
imaginary (euclidean) time, it becomes two surfaces, future and past horizon in real 
time. In this section, we show that the singularity appears only in the past horizon 
even for 7 ^ 0, if we take the Eddington-Finkelstein coordinates. Here, we focus 
on the near horizon region and discuss about the regularity of the gauge held at the 
horizon. 

In the near horizon region, the metric of the non-extremal black holes is univer¬ 
sally given by the Rindler space; 

ds^ = —rdt^ H-h {dx'Y , (E.l) 

r 

where r = 0 is the horizon of the black hole. Only outside of the horizon is covered 
by this coordinates. In order to move to the Eddington-Finkelstein coordinates, we 
dehne null coordinates as 

t± = t ±\ogr (E.2) 

and then, the metric is expressed as 

ds^ = —rdt\ ± 2drdt± + {dx^Y ■ (E-3) 

The ingoing (outgoing) Eddington-Finkelstein coordinates (with t_|_ (t_)) also cover 
inside of the future (past) horizon. The Kruskal coordinate is dehned by 

and then, the metric becomes 

ds^ = dx+dx- + (dx*)^ . (E.5) 
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This covers all region, and a;+ = 0 and a;_ = 0 are the past and future horizon, 
respectively. 

If the gauge held has non-vanishing and regular A,, at the horizon, it is singular 

there. It can be seen as follows. In the Kruskal coordinates, the gauge held becomes 


A — A^dt 




(E.6) 


Therefore, the gauge held is singular at the future and past horizon. In the Eddington- 
Finkelstein coordinates, it is expressed as 

T -) . (E.7) 

r J 


A = A, 


However, if we take the ingoing Eddington-Finkelstein coordinates, and if Ar is 
not singular at the horizon, the gauge held is singular only at the past horizon and is 
regular at the future horizon. In the Kruskal coordinates, the gauge held is expressed 
as 

A = A^dt+ + A^dr = A^—— + Ar (x+dx_ + x^dxjA . (E. 8 ) 

x+ 

Therefore the gauge held is singular at the past horizon = 0 but regular at the 
future horizon a;_ = 0 . 


F. First order solution for general 2 ; 


The correction terms can be calculated in a similar fashion to the z = 2 case. We 
dehne h^u, O/i and (p as 

gu = r^^{-f + hu) , (F.l) 

1 

9tr — 2 ’ (E-2) 

gti = r'^hti , (F.3) 

gij ^ hxx^ij T hjj) , (F.4) 

A = [a{x) (r® —''"o(^)) — Ai{x)v'^{x)~\ dt — a{x)r‘^dr + Ai{x)dx'^ + atdt + aidx"^ , (F.5) 

and 

^ = 4>-^ log(/ir"®) . (F.6) 

In the n*(0) = 0 gauge, the constrains are expressed as 


p{x) = a(a;), (F.7) 

dta = adiv'^ , (F.8) 

dtTo = —rodiv'^ , (F.9) 

dtAi + AjdiV^ = " 2 !^ ^| rg+^a(9jro . (F.IO) 
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The first order solution for the sound modes is 


htt = 2 \/ 6 - 


z + 3 




2^-(;.-3)^ hi' 


!£!1 .(•) 


+ ^ ^ {z + sy 


ro^+n 


/ a 


+ r-"-" 1 - 


h2(r) — r ^ ^ 2^: — ( 2 : — 3) 


dr h 2 (r) 


V3Z“|-3 o yZ ~\- 0 ) Oj 


dr h 2 {r) 


(F.ll) 

(F.12) 


at = a® — 3 ^ ^ ah^i^ -^+ \ — - — ar^~^^c 

z — 1 z + 3 \ z — 1 

- ar^div'' 

3 

, . ,<«) _ y^,.-.-3,40, _ 


+ 3 ^ [ dr h 2 {r) 


z-1 


(F.13) 


3 1 


r ^ ( rh 3 (r) — 2(;2 — 1) / dr h 2 {r) 


(F.14) 


where 0®, hf'\ a® and are integration constants. The function ^2 is the 
solution of the following differential equation; 


0 = (r^+^ - rZ^) Kir) - r [zr^+^ + {2z + 3)r^+^] Kir) 

-iz + 2) [i2z + 7)r^+^ + iz + 2)rK^] h 2 (r) 


(F.15) 


and is given by 


h2ir) = 2 ^_. 

+ C'2rK^+i-VF+3)(9-+i9))2Fi 


(F.16) 


where Ci and (72 are integration constants, and 


Q:± = 


z + 3 


The hrst order solution for the vector modes is 


hti — 


hi { 4 ") _ 

r6-^ \ a 


(F.17) 


(F.18) 


= ( 2 ( 2 ; — — ( 2 ; — 5)rQ+^) Jdraiir), (F.19) 
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where and a. 


,(o) 


are integration constants. The function Oi is given by 


a, (r) =__ 

^ ^ (r"+3 - rl+^)[2{z - l)r"+3 - {z - 5)r^+3]2 

a2(r) = C^-{z + 3)ar"-^ (r"+3 _ ^-+3) /.W 

7 -I- .^^2+2 

+ a (10(;^ - l)r^+3 _ _ 5)^^ _ 2)^.+3) 


(F.20) 


(F.21) 


and 02 (r) is expanded around r = tq as 


/ N r- 1 z{z + ?,farl^dirQ 
a2(r) = (Fa +- 

In order for the solution to be non-singular at r = tq, we have to take 

z{z + 3)^arg^(9iro 


(F.22) 


C'a = 


2 ( 2 - 1 ) 


(F.23) 


and then, 02 (r) becomes 

a2(^) = —{z + 3)ar^~^ ^^^+3 _ h® 

r -I- 

-\— - --^^a ( 10 ( 2 : — l)r^’''^rQ + z(z + 3)r^rn — (z — b){z — 2)r^~^^') dir^ . 

2(;2 — 1) ^ 

(F.24) 


The hrst order solution for tensor modes is 

r‘^dr 


f~tjA - (T j ■ 


+ ^4 


dr 


/ ^ 2+3 _ ^2+3 ' J ^(^^ 2+3 _ ^ 2 + 3 ^ 


(F,25) 


The regularity at tq implies 6*4 = Tq. Then, we obtain 


r(r.+3 - ryl) ' 


(F.26) 


G. Counter terms for general 2 ; 


In order to obtain regular stress-energy tensor, we introduce the counter terms; 


= 


IGttG 


d'^XyJ—''^ 


-{A + 2z) + C+ 

-*-1 


(G.l) 
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Then, the renormalized stress-energy tensor is obtained as 


rp\j _ 

0 — 


T\ = 


T°- = 


_ 


1 


SttG 
1 1 
SnG 

1 z-1 

SnG ar^+3 
1 1 


G z + 2 
¥ ^ 


^,+3 _ f 


z + 3 


a 

z{z + 3) 


z-1 


(G.2) 


rQ^diVo — - - -v'^v^Aj + ^rlcTijV^^ + Oiv' 


^{z-l) 


a 


Ai + C>(r-(^+^)) 


i(l + C)rS*%, - h„Vy + ^«'A, 




^ SttG r"+3 

Since for general ; 2 , the volume form on the boundary behaves as 




‘f 


,z+3 


(G.3) 

(G.4) 

(G.5) 

(G.6) 


the above 0{r ^ terms are regular. For G = z — 1, the energy density and the 
pressure becomes 


T = — 


P = 


y^^ + 3 


(G.7) 


levrG ° ’ ^ 167rG' 

and they satishes 

zS = {d-l)P. (G.8) 

The constant G in the counter terms can be hxed by the regularity of the dual 
operator to the dilaton 0. By introducing the counter terms, the expectation value 
of the operator becomes 


Os = 




IQtiG 

+ . (G.9) 

In order for the above expression to be regular, we have to take G = z — 1, and then, 
we obtain 


T°n = 


TG = 


1 1 
SttG 
1 1 


.3 .+3 _ £— \ ^ 
2 ° a 


A^ + , 


(G.IO) 


z + 3 


° SttG 

+ C>(r-(^+4) 


z{z + 3) 




4{z - 1) 


f”, = ht G , + Cl(r-l*+‘‘l) , 


8ttG ar^+3 

. — — _3_ ( _ -r^a - + ^ ~ 


T\ = 


and 


o^ = 


IGttG 




(G.ll) 

(G.12) 

(G.13) 

(G.14) 
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